Topological Susceptibility and Zero Mode Size in Lattice QCD 
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We use the overlap formalism to define a topological index on the lattice. We study the spectral flow of 
the hermitian Wilson-Dirac operator and identify zero crossings with topological objects. We determine the 
topological susceptibility and zero mode size distribution, and we comment on the stability of our results. 



1. INTRODUCTION 

In this proceedings, we summarize our work on 
the determination of the topological susceptibil- 
ity and zero mode size distribution of SU (2) and 
SU(3) gauge fields using the overlap formalism. 
Most of this work has appeared in print jlj . We 
will focus here on new results which are exten- 
sions of the published work, namely the topolog- 
ical susceptibility for SU(2) pure gauge theory, 
the zero mode size distribution for both SU(2) 
and SU(3), and the density of zero eigenvalues 
p(0) for SU(2) and SU(3). 

A more detailed discussion of the role of 
the zero eigenvalue density can be found in 
the proceedings of the plenary talk given by 
Narayanan ||. 

2. SPECTRUM OF THE HERMITIAN 
WILSON DIRAC OPERATOR 

The overlap formalism for constructing a chi- 
ral gauge theory on the lattice ||] provides a 
natural definition of the index, /, of the associ- 
ated chiral Dirac operator. The index is equal 
to half the difference of negative and positive 
eigenvalues of the hermitian Wilson-Dirac opera- 
tor Hl(to) = 75 W(— m) where W(m) is the usual 
Wilson-Dirac operator (we will use negative sign 
for the mass term throughout). 

On the lattice, because of the additive mass 
renormalization, the crossings of zero occur at 
positive m and spread out in m. It is easy to 

'Speaker at the conference. This research was sup- 
ported by DOE contracts DE-FG05-85ER250000 and DE- 
FG05-96ER40979. Computations were performed on the 
QCDSP, CM-2 and the workstation cluster at SCRI. 



see, that no eigenvalues of Hl (to) can be zero for 
m < 0. Since in the free case the first doublers 
become massless at m = 2 we restrict ourselves to 
m < 2. A simple way to compute the index / is to 
compute the lowest eigenvalues of Hl(to) at some 
suitably small to before any crossings of zero oc- 
curred. Then to is slowly varied and the number 
and direction of zero crossings are tracked. The 
net number at some mt is the index of the overlap 
chiral Dirac operator. 

We have applied this procedure to compute the 
index, which we take as the definition of topolog- 
ical charge on the lattice, of various gauge field 
ensembles with gauge group SU(3) including pure 
gauge, Symanzik improved pure gauge, and two 
dynamical flavor generated gauge fields Q), and 
recently pure gauge SU(2). We found that the 
zero crossings start occurring at some, ensemble 
dependent, toi > and continue occurring for 
all to in mi < m < 2 in sufficiently large lattice 
volumes. Hence, the spectral gap is closed in the 
entire range. 

To further investigate the zero level cross- 
ings, we compute the density of zero eigenval- 
ues p(0;to) for Hi(m) by fitting the integrated 
density p(X') d\' to a line through the origin 
for some cutoff A. By studying the scaling of 
p(0; m) for quenched SU(3), we conclude that the 
density of zero eigenvalues falls exponentially in 
the inverse lattice spacing and only vanishes in 
the continuum limit (see Figure 6 of Narayanan's 
talk Q). We will observe latter that these zeros 
are due to small localized modes. More discussion 
can be found in the talk of Narayanan B . 
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Figure 1. Detailed study for = 6.0, 16 3 x 32. 
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3. TOPOLOGY AND SMALL ZERO 
MODES 

Using the index of the chiral Dirac operator 
as our (fermionic) definition of the topological 
charge of the gauge field background, we obtain 
the topological susceptibility x as a function of 
m t . All of the gauge ensembles studied have the 
general characteristic that the susceptibility rises 
sharply in the region where p(0; m) is peaked, and 
then it essentially flattens out in the region where 
p(0; m) is small [0. 

We define a size of the eigenvector associated 
with the level crossing zero mode 



Pz 



/(t) = £tr(0t(M)^(M)) 

^ /max 



motivated by the t'Hooft zero mode p1/2(t 2 + 
Pz) 3 ^ 2 where cf>k are the eigenvectors of 
Another definition based on the second moment 
of f(t) was used in Ref. JO. We should empha- 
size that we look only at the sizes of eigenmodes 
that cross, and only close to the crossing point. 
Only then can we expect to get a good estimate 
of the localization size inspired by the 't Hooft 
zero mode. 

We show in Figure |l| a detailed study of the 
(3 — 6.0, 16 3 x 32 pure gauge ensemble. On the 
first line is shown the density of zero eigenvalues 
p(0; m) and the number of crossings in each mass 
bin. Since there are a nonzero number of cross- 
ings, we see that p(0; m) does indeed measure zero 
eigenvalues, and not just small eigenvalues near 
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Figure 2. Zero mode size distribution. 
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Table 1 

Topological susceptibility and parameters. 



scale. Our results are in rough agreement with 
other groups P~|6[. 

In Figure Epwe plot the zero mode size distri- 
bution for several /3's for SU(3) and SU(2) which 
is of interest for the instanton liquid model. We 
find the distributions always peaks for sizes about 
1-2 lattice spacings. For 5/7(3), we see small 
secondary peaks at ~ 0.65fm for (3 = 6.0 and 
5.7. However, there is only a shoulder for 5.85 at 
0.3fm. For SU(2) which are all the same lattice 
size, we find the secondary peaks shift to smaller 
p z for increasing (3. These secondary peaks do 
not occur at any general physical radius and are 
consistent with a finite volume effect. 



zero. We also see that p(0;m) rises sharply in 
to, then falls to a nonzero value where there is a 
small number of zero level crossings. 

On the second line of Figure [l], we show the size 
of the zero modes p z (m). The modes are large 
near mi where p(0; to) is large, then p z drops to 
about 1 or 2 lattice spacings up to m = 2. We 
see that the corresponding x rises sharply when 
p z is large for to near mi, then is quite stable 
when p z is small. This result show that while the 
index, /, of the field is to dependent, \ ( a physical 
quantity) is independent of the contribution from 
the small modes for m > 1. 

To further clarify the relative contribution of 
the zero modes, in the last line of Figure the 
zero mode size distribution is shown which peaks 
for p z < 2. In the adjacent graph, x, here defined 
by the contribution of zero modes of size p z and 
larger, is stable when p z < 2. Hence, the small 
modes do not affect the estimate of x- Our esti- 
mates of x are shown in Table [l] where we use the 
string tension value s/a — 440(MeV) to set the 
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